clear;
clc;
global mu;
global lambda;
global b;
global c;
global a;
global eta;
global K;
theta=0.1:0.1:pi./(2.*a);
global phita;
global stepsize1;
stepsize1=2.*pi./(pi./(2.*a)/0.1);
phita=stepsize1:stepsize1:2.*pi;
a=1;
b=0;
c=0;
ed=(1+(1./4).*b.*(3.*cos(2.*a.*theta)+1)+3.*c.*sin(a.*theta).^2.*cos(2.*phita)).*a.*sin(a.*theta)./(2.*pi.*sin(theta));
theta
ed
phita
[X,Y,Z]=sph2cart(theta,phita,ed);
% surf(theta,phita,ed);
surf(X,Y,Z);
